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perform nonlinear constrained optimization of engineering 
design problems. The program is developed especially for 
use Oh wicrocomputers and is called Microcomputer Software 
for Constrained Optimization Problems (MSCO¥). It ya 
accept a nonlinear objective function and up tc ae 
inequality constraint functions and up to 20 Founded design 
variables. 

MSCOP employs the method of feasible directions. 
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ment, the MSCOP was implemented on an IBM 3933 uSing stan- 
dard Easic language, Waterloo BASIC Versione Tt aes 
directly transportable to a variety of microcomputers. 

Typical applicaticns of MSCOP program are in the design 
of machine components and Simple beam and truss structures. 
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I. INTRODUCTION 


A. PURPCSE 


This thesis describes the development of a micrecomrfuter 
oriented program called MSCOP (Microcomputer Software for 
Constrained Optimization Problems) for constrained optimiza- 
tion of engineering design protlems. Problems which can he 
solved ty the MSCCF are nonlinear programming frotlems 
arising in several areas of machine and structural design, 
such as the minimum weight design of structures subject to 
stress and displacement constraints [Ref. 1}. 

in recent vears, several powerful general purpose ovti- 
mization programs have become available for engineering 
design problems, @.g., COPES/CONUD GY) Reraap and ADSoa 
PRe we. oS ie These prcegrams can hanile a wide range of design 
problems and contain a variety of solution technicués. 
Also, several programs are available that include optiniza- 
tion in an integrated analysis / design code, e.g., ACCESS, 
ASO? ,- EATS PARS SAVES, SrAn, STARS and TSO [Ref. 4]. Alar 
of the akove optimization programs are written in FORTRAN, 
and are Euilt for use on a mainframe computer. Their use can 
be cumbersome, especially for the occaSional user. Since 
Many engineers are now using microcomputers, there is a need 
to develop an optimization program contained in a microcene 
puter software package for use on microcomputers. This 
thesis fills that need by developing a compact program 
written in a standard BASIC language suitable for a wide 


range of microcomputers. 
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B. I EPLESENTATION 


The nature of an optimization program depends on the 
computer and programming method available. Tie SCOR sorts 
ware is designed for use on a microcomputer. However, <0Or 
the speed of development and testing, MSCOP was developed on 
mye 6EM 3033 computer at the ¥. R. Church Computer Center in 
Naval Pcstgraduate School, and was written in WPASIC 
(Waterloo BaSic) Version 2.0. 

To make sure that the program is easily portable toa 
micrccomruter, only standard BASIC commands and functions 
@ee used. For example, FOR JT = 1 TO IDB ... NEXT I, GOSUB 
etc., were used. The commands and functions not availatle 
in all variations of EASIC are avoided, for example, TRN(A), 
MAT(A), etc. 

MSCOP provides design engineers with a convenient tcol 
for optimization of engineering design problems with up to 
20 bounded design variables andas many as 50 ineguality 


constraints. 


C. GENEKAL OPTIMIZATICN MODEL 


The general optigization problem to be solved is of the 


Zatee bind athe Set of deSign variables X that will 


Minimize Eee) (Veet) 
sup ject to eG) ee J = 1,--6, (lee) 
J — 
al u 
Mew. Sk Aen earl pratoxcigne 1 (Teo) 
a ren 3h 


where X is referred to as the vector of design variatles. 

mi.) iS the objective function which is to be minimized. 
L 

eye are inequality constraint functions, and X; and i 


are lower and upper rounds, respectively, on the design 
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Variables. Although these LEounds or "side constraints" 
could be included in the inequality constraint set given by 
EFqa(1.2), i1t 1S convenient to treat them separately because 
of their svecial structure. The objective function and 
constraint functions tay be nonlinear, explicit or impligeaee 
ee However, they must be continuous and should have 
continuous first derivatives. 

In general engineering optimization problems, the okjec- 
tive to kre minimized is usually the weight or volume ofa 
structure being designed while the constraints gives limits 
cn compressive stress, tensile stress, Euler buckling, 
displacement, frequencies (eigenvalues), etc. f Ref. (ame 
p-264]. Equality constraints are not included because their 
inclusicn complicates the solution tecanigues and because in 
engineering Situations, equality constraints are rare. 

Most optimization algorithms require that an initial 
value of design variables X° be svecified. Beginning from 
these starting values, the deSign is iteratively improved. 
The iterative procedure is given by 


+1 
ie x + ax s? (1.35 


where gq is the iteration number, S 1S a search direction 
vector in the design space, and a* iS a Scalar parameter 
which defines the amcunt of change in X. At iteration gq, it 
is desirable to determine a direction S which will reduce 
the objective functicn (usable direction) without violating 
the censtraints (feasible direction). After determining the 
Search direction, the design variables, X, are updated Ly Eq 
(1.4) so that the minimum objective value is found in this 
GIECCt1Onus di he tego. 

Thus, it is seen that nonlinear optimiZation algonicime 
for the general optimization proklem based on Fq(1.4) can be 
separated into two parts, determination of search direction 


and determination of scalar parameter a*. 
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Dee OnGANTZATION OF THiS THESIS 


This chapter has stated the purpose of the thesis and 
has put the general ccncept of engineering optimization into 
a preliminary perspective. Chapter 2 will describe the 


essential aspects of the optimization algorithm used in 


MSCOP such as finding a search direction, the one- 
dimensional search and convergence criteria. Chapter 3 
describes program usage. In chapter 4, there are three 


examples which are sclved by the MSCOP. Summary and conclu- 
sions are given in chapter 5. The program is listed ir the 


appendix. 


re 


IIT. CPTIMNIZATION ALGORITO 


= ee SS ee ee ce ee a SS SS ee ee 


A. INTRODUCTION 


There are many optimization algorithms for constrained 
nonlinear problems such aS generalized reduced gradient 
method, feasible direction method, penalty function methods, 
Augmented Lagrangian multiplier method, and seguential 
linear programming. The feasible direction nethod is chosen 
for development in this thesis for three main reasons. 
First it progresses rapidly to a near optimum design. 
Second it only reguires gradients of objective and 
constraint functions that are active at any given poirt in 
the optimization process [Ref. 7]. Third, because it Maas 
tains a feasible design, engineer cannot fail to meet safety 
reguirements as defined by the contraints. However, the 
method does have several disadvantages in that it is fprone 
to "zig-zag" between constraint boundaries and that it is 
usually does not achieve a precise optimun. This method 
solves the nonlinear programming problem by moving froma 
feasitle point (can be initially infeasible) to another 
feasible point with an improved value of the objective 
value. 

The following strategy is typical of feasible direction 
method : Assuming that an initial feasible point Xoo 
known, first find a usable-feasible direction S. The algo= 
rithm for this is similar to linear programming and Compaes 
mentary pivoting algorithms. Having found the search 
direction, a move is made in this direction to update the X 
VeECtor according to Eqa(i.u)a. The scalar a*®* is found Dye 
one-dimensicnal search to reduce the objective function as 


much as possible subject to constraints. That is |i 
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ete. 
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PRINT RESULTS 


CALI = CALI+1 


no 


— 
<CA11 = as 
oO 
10% RENUCE THE 
DESIGK VARIABLES 


OCI Ne Rease THE 
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Figure 2.1 Algorithm for the Feasible Direction Method. 


Peerats) Subject to G{xXta*S) < 0. It is assumed that the 


initial design X° is feasible, but if it is not, a search 
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direction is found which will direct the design Toma 
feasible region. After updating the X° vector, the conver- 
gence test must be performed in the iterative algorithm. A 
convergence criteria used in this is implementation are 
described in section f. The general algorithm used in “SCOP 


iS given in Figure 2a) 


Be SEARCH DIRECTION 


In the feasible direction algorithm, a usable - feasible 
search direction 5 is found which will reduce the objeg. 
tive function without violating any constraints for some 
finite move. It is assumed that at any point in the design 
Space (at any X) the value of the objective and constraint 
functions as well as the gradients of these functions with 
respect to the design variables can be calculated. Since 
these gradients cannot usually be calculated analytically, 


the finite difference method Eq(2.1) is used in MSCOP. 


ar (ZX) F(X+ €e.) - F(X) 
i ape ee a (2am 


aX. E 
a 


where e is the ith unit vector 
a: 


&is a small scalar. 
In MSCOP, € is 0.1% of the ith design variable 


In the feasible direction aigorithm, there are usually 
one Or mere “active" constraints. A constraint G(X) 
“active" (at K 22g (x) ea or As shown in Figure 2.17, lima 
constraints are active the standard steepest descent direc- 


CLOT SS) see Vales sca. 
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1. USable-Feasible Direction 


a a A an SE cm: cm en a a ea a RS 


F(X) =const 
Xa(X°) 


epee 


sector 





/) Usable-feasible 
47. sector 





Figure 2.2 Usable-Feasible Direction. 


Assume there are NAC active constraints at X. The direction 


Sis “usable" if it reduces the objective function, i.e., 
Whi Seceae (2. 2) 


Similarly the directicn is feasible if for a small movement 


in this direction, no constraint will be violated, i.e., 
Wer Sie< 0 (2. 3) 


This 1s shown geometrically in Figure 2.2 


by 


2. Active Constraints 


It is necessary to determine if a oconstraint is 


active or violated in the feasible direction algorithn. A 
constraint G(X) ¢ 0 is "active" at X9° if G(X°) =0. in coaee 
to avoid the Zigzagging effect between one Or more 


constraint boundaries, a tolerance band about zero is used 
for determining whetker or not a constraint 1S active. From 
the engineering point of view, a constraint G(X) <¢ 90 is 
active near the boundary G(X) = 0 whenever ACC € G(X) £ VCC. 
AcC is the active constraint criterion and VCC iS Gime 
violated constraint criterion Win @%seor- Assuming the 
feasikle constraints are normalized so that G(X) ranges 
between -1 and 0 for reasonable values of X, the constraint 
G(X) < 0 1S considered activo: yo > -0.1. The 
constraint 1S considered to be violated if G(X) > 0.C0ae8 
This is an algorithmic trick which improves efficiency ani 
reliability of the alcorithn. However, Since in the one - 
dimenSicnal search, all interpolations for constraint G(X) 
are done for zeros of a linear or quadratic approx imaticn te 
G(X) in erder to find a*, at the optimum the value of active 
constraints are very near zero, but may be as large as 0.004 
fRef. 6]. From an engineering point of view, a 9.4 fF 


constraint violation is considered to be acceptable. 


Zoutendijk [Ref. 8] has shown that a uSabDtouee 


feasikle direction S may be found as follows 


Maximize 8 (2.55 
Sup mec t. tow 


GF(X)-S +@< 0 (25s 


TAG jaort 9 3 | < 0 jeg (2. 6) 
S 
S bounded Za 


Where scalar Bis a measure of ficmesdtersraction of the 
usability and feasibility requirements. The scalar 6; in 
Eq (2.6) is referred to as the "push-off" factor which effec- 


tively pushes the search direction away from the active 


es ae Ne cen 





Pigure 2.3 Push-Off Factor and Bounding of the S-Vector. 


Senstraints. Pieces O)ip te che = PUSh-oLk Lactor 1S Zero, 
the search direction is tangent to the active constraints, 
and if it is infinite, then the search direction is tangent 


zo. the objective function. It has been found that a 
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Fush-off factor is defined as follows gives good results 
(Ret. 255. Pp. i677 |e 


bo 


@ =| 1-—— O. (2. 8) 


where @ =1. 


To avoid an unbounded soltttion when seeking a usable 
- feasible direction it 1s necessary to impose bounds on the 
search direction S. Cne method of imposing bounds on séarch 
Girection iS to impose bounds on the components of S-vector 


C25. Orm x 


= ncaa (29h 
This choice of bounding the S-vector actually biases the 
search direction. This is undesirable since we wish to use 
the push-off factors aS our means of controlling the 
search direction. A method which avoids this bias in search 
direction is the circle as shown Figure 2.3 . The norm here 
1s 

s°3 <1 (2.925% 


— 
7 


4. Simple Simplex-like Method for Search Direction 


— = -— —_— —_— —_— — ee eee coe cee — a oe oe —_— — a ee ce — a ee 


Vanderplaats [Ref. 5: yes ee lS] provides the 
matrix formulation which solves the above sub-optimization 


rroblem by using the Zoutendijk method. 


Maxinpize Bay (2. 1G 
SUubjeectmto: ; 


A-y <0 (2. 108 
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I< 

I< 

LA 
= 


Where 


lM {WH 
=e 


Dn 
7] 
i el 


AB 

eee 8 
9G. (x 

aA > | de 8. 


/ : 
7G (xX) 
oe eee 


a 
VE a(x), 1 


es 


Pk 1B) 


ity 
il 
= GC) ee © co) ©} 


C22 14) 


and where j is the number of active constraints (NAC) 


Mien the solution to Eq(2.10) through (2.12) is found, S$ may 


be normalized to some value other than unity, but the form 


of tne normalization is the same. Meso lution to the above 


problem may be obtained by 


derived from the Kuhn-Tucker 


Where 


solving the following system 


COnult TONS, fOr that probienr 


(2. 15) 


1a 


0 ee (22S) 


Ze 


TT 


eae (2. 17) 
iF Cememey Ma easy (2. 1s 
GS S018 (2. 125 


ALove system can be solved using a comnlimentary pivot algo- 
ieaat alee Choose an initial basic solution to Ea({(2.15)” dome 
ke 


Vo= 3c; Wi = 0 (2. 265 


where v is the set of basic variables and u is’ the set of 
nonbasic variables. If all v; > 0, Eq(2.16) is alsousamaee 
fied and problem is solved. If some v; < 0, the solutagg 


procedure is as follews : 


Let BE; be the diagcnal element of the i-th nonbasic vam 
able. 

1. Given the condition that some c is less then zero, 
we find max (c;/Bii) which is the incoming row to the 
basis. 

26 whe J ne Oma: 1G column is changed to a kasic 
colunn, the tatleau is updated by a standard simplex 
PIVvOle. On Ba. 

3. Untal all ¢;> 0, fepeat steps tieauce. 

Po. Wen ease oc, the iteration 1S complete. The 
value of u is now the desired solution. 

a2, Syeising ee p-a‘u, we get the usable-feasible 


search directicn S which is first NDV components of 


Ye 


The method of feasible directions assumes that we 
begin with a feasible design and feasibility 1S maintained 


throughout the optimization process. If the initial design 


a 


is infeasikle, then a Search direction pointing toward the 
feasible region can ke found by a simple mnodificaticn to 
direction finding prctlen. 

A design Situation can exist in which the violated 
constraints are strcngly dependent on part of the design 
“variables, while the objective function iS primariiy depen- 
dent on the other design variables. This suggests a method 
for finding a search direction which will’ simultaneously 
Minimize the objective while overcoming the constraint 
violations. These considerations lead to the follcwing 
statement of the direction finding problem [Ref. 5 : 
epee? I-72) : 


Maximize eer )e* ee BBR (Zee 
Subject to ; 


yG(X)-S+ 868 <0 j€g (2. 22) 
j 


S°S < 1 (2. 23) 


where J is the set cf active and violated constraints, and 
where the scalar @ in Fq(2.21) is a weighting factor deter- 
Mining the relative importance of the objective and the 
constraints. Usually a value cf $& > 10000 wili ensure that 
the resulting S-vector will point toward the feasible 
meeron. Incorporating Eqg(2.21) and Eq(2.22) into the direc- 
tion finding algorithm requires only that we modify the 


p-vector given in Eq (z.24) and the A-matrix of Eq(2.25). 


- WF (x) 
p= (22028) 
Ey 


DES: 


Gace 
A = ¥ 2' 9» (2.25) 
p 
VG (Xx), 
Y O 
Garreaeo (2 26% 


We use the Simple simplex-like metnod to find the 


search direction toward the feasible region. 


Ce CNE-DIBSENSIONAL SEARCH 
1. No Violated Ccnstraints 


TF no constraints are violated, we find the largest 
a*¥ in Eg (1.4) from all possible values that will minimize 
the objective on S without violating any constraints, active 


Or Inactive. 


The procedure in MSCOF is as foliows ; 

1. Let a0, al, a2, a3 be the scalar in 3q(1.4%) Ge@pnee 
sponding to peints X09) x19 27k 3 eee 

20 “ath =)0 at Gay ch pean teed: 

3-5 DNROEdery SvOe dct al, we can calculate the al to 
reduce the objective by at most 10% or to change each 
of the design variable xX by at most 10%. 

4. Update the design variables to X1 uSing Eq9(1.4). 

5. Evaluate the cbkjective for Xl, and check the feasi- 
bility. If one or more constraints 1S violated, jegea 
al is reduced to al/2, and we go to step 4. 

6. In order to estimate a2, we can use the quadratic 


approximation with 2 points X,) ks eandee eye 


feaeicdaterthe = GesiGgn varlables to £2 by Eq(1. 4) and 
check the side constraints. 

8. Evaluate the objective and constraints. 

9. Now having 3 a's, and values of objectives and 
ccenstraints for design variables XQ, Mle, “X25 are 
LOM cOmnye Using S—-polnt Quadratic approximation, a 
value of a3 is found. 

10. Update the new optimal point in search directicn by 
Eg (1.4). 

11. Evaluate the okjective and constraints. 

12. Now choose last 3 values, al, a2, a3and find a new 
qoelemnges-POINntS OUadratic approximation 

13. Choose the a* among the 5 points which corresfonds to 
the minimum objective function value with no-viclated 


Gensttralnts. 


Ponce Omeiiore Comsttaints Violated 


If one or more constraints are initially violated, a 
modified usable-feasible direction is found. It is then 
necessary to find the scalar a* in Eq(1.4) which will minri- 
mize the maximum ccnstraint violation, using the nost 


violated constraint j, a good initial estimate for a* is 


~G (K) 
ax = ———____ ones 


Since the gradients of the violated constraints are 
known, the scalar which is reguired to obtain a feasitle 
design with respect to violated constraint in the search 
BeecCetION, 1S given to a first approximation by Fq(2.27). 

The more detail procedure in 4SCOP is as follow ; 

1. Choose the most violated constraint j. 
2. Calculate a* for violated constraint 5 using 
MON 2s27,). 
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3. Update the design variables for a* and check the side 
COns tla dn esi. 

4“. Tf one or more violated constraints still exist, then 
calculate the derivative of objective, violated and 
active constraints and find anew search direction 
and then go .te Step ai. Otherwise proceed with the 


optimization in the normal etaciwiornm 


BD. CCNVERGENCE CRITERIA 


A desired property of an algorithm for solving da none 
near proklem is that it should generate a sequence of points 
converging to a global optinal @eperm ce In many cases, 
however, we may have to be satisfieil with less faverable 
outcomes. In fact, as a result of non-convexity, piceien 
size, and other difficulties, we may stop the iterative 
procedure if a point kelongs to a described set, which is 
defined in MSCOP as fcllows ; 


ond 
a 
tO 
li 
nm 
{>< 


1 1 eo sy eae 


NO 
tO 
iN 


LU Ee te eerie race ee 


In MSCOP, the algorithm is terminated if a point X is 
reached such that XEQ, (1) Oe &,1s 0.091) amd ee is 
approximatly 0.001. Since in engineering design problems it 
is not necessary to find solutions with more than three 


SIGniireene Gidi ts. 


26 


Pete "HSCOP USAGE 


A. INTRODUCTION 


Sance thas | MoCOr 2S written in WATERLOO BASIC Version 
eel tt) 1S yery convenient to use. The user must first 
formulate the design problem with the classical machine 
Cesign criteria. Given the formulation of the design 
problem as a nonlinear progran, the user then enters the 
problem aS a part of a BASIC progran. The user defines the 


Sepective. f£uUnetion and constraint functions using EASIC 


statements. Other parameters are input as data : the number 
of design variables NDV, the number of ineguality 
@enstraints NIOQC, variable bcunds an initial design value 


or a PrIakt Control number. 


B. PRCBLEM FORMULATICN 


Generally, the experienced design engineer will be able 
to choose the apfropriate objective £fOEmeeenp Limi Zation 
depending on the requirements of the particular application. 
The physical phenomena of Significance should first be 
summarized for the device to be designed. The appropriate 
objective can then be selected and constraints can be 
imposed cn the remaining phenomena to assure an acceptatle 
design from all standpoints. However, the initial fornmula- 
tion for the optimization problem should not be more compli- 
cated then necessary and this often requires the making of 
some simplifying assumptions. [Ref. 9]. 

After completing the formulation of the design ctrctlen, 
the design engineer should be able to answer the following 
questions 


1. What are the design variables ? 
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2. What is the objective function ? 
3. What are the inequality constraints ? 
4, What are the bounds on the variables ? 

The engineer is then ready to input the program to the 
MSCOFP. However, additional study and preparaticn of the 
problem may be useful. In particular, redundant constraints 
snould be avoided if possible. MSCOP will operate with 
redundant constraints but it will operate faster without 
then. Selection of an initial deSign point from which 
Start this program is important, Since it affects peneoe. 
ance and running time. The user should use any availatle 
information which gives a good initial approximtatvom is 
Side constraints exist, the user must be sure the initial 
values cf the design variables do not violate the side 
CONS traits. This frogram will automatically handle an 
initial design point which is infeasible with respect to the 
G(X) < 0 constraints. However, if the initial point does not 
violate these constraints, convergence will likely be more 


Ea pale. 


Co PRCBIEN SWary 


Froblem entry is accomplished by editing the tain 
program directly. As an example, consider the following 
Simple NLP with two design variables, and three constraint 
Linctaons. 


Z 2 
Minimize F(X) = X + 3 KK. 9-472 ee ee 
1 ieee 2 1 2 


SUbHIECe EOS 


1 1 
— +—-2<0 
xX X - 
1 2 
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With the MSCOP loaded into memory and listed on the CRT, 
modifications are made on the program lines as’ follows to 


input this example : 


Line 109 
Just after the word "data", three integers are added, 
separated by a comma. The first number is NDV which is 


the number of deSign variables, the second is NI¢C which 
is the number of inequality constraints, and the third is 
Peiiewrch is DEint Control number (0 : only firai 
results, 1; given data and final results, 2 ; given data 


and iterative subcrptimal results) 


for example : 
WOON data s253, 2 


Lines 291-220 


Each line here corresponds to a separate design variable, 
beginning with X(1) and continuing in order to input 
X (NDV). On each line, three values are separated by 
commas. After the word "data", these values are the 
initial values of the design variable, the lower bound on 
the variable and the upper bound on the variable. iE eno 


bound 1s to be specified, the entry is filled by "no". 


For the sample protlem, the input is : 
ZO wedata Se, 0. 1,10 
PUCo data. .3., 0. 1, no 


Zo 


Lines 400 - 450 


These lines are available for defining the objective 
LUNG tren. The ckjective function must be defined in 


terms of subscripted design variables X(1), X{(2), etc. 
For the sample proklem, the vmyiiteoee 
4OO f£n_f = x (1) ¥**2+x(1)*xX(2) +2. *x (2) **2-x (1) -x (2) +1. 
tanes” 500-6590 


These lines are available for defining the inequality 
constraint functicns, which must be expressed using the 


format : 


601 1f 2 = K ther (nige— Cece 
1 at 


For the sample protlem, the input is ;: 


00601 if 1 = 1 then fn_g = rea eas Se 
O0602 26 1°= 2 tl ete G ae SX (242 
G90603 1f 1 = 3 then Eng = x (1) €¥2+x% (1) -x (eee 


If there are many constant values in the constraint func- 
tions, the user may input data for these functions on 


lines 501-600 in order to simpiify their statements. 
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iV. EXAMPLE PROBLESS 


A. DESIGN OF CANTILEVERED BEAM 
(Ui comteeatiit levered Ream 


Assume a cantilevered beam as shown in Figure 4.1 


must be designed. The objective is to find the rininun 


20000 Psi 


es] O} 
T ii 


50 220m bs 1 
Ve ole O ire m 


P = 10000 lbs 





ae 


Figure 4,1 Design of a Uniform Cantilevered Bean. 


volume of material which will support the load P. 
The design variables are the width B and height H in 
the team. The design task is as follows : Find B and H to 


tinimize volume V =B Fl (4. 1) 
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we wish to design the beam subject to limit on bending 
stress, shear stress, deflection and geometric conditions. 


The bending stress in the beam must not exceed 20,000 psi. 


Me 6 P 1 
G = = (4. 2) 
b 2 - 


Bore 





The shear stress must not exceed 19,000 psi. 


3 P ee 
0, = = ——— < 10,900 (4. 3) 
2 A 2B H = 





and the deflection under the load must not exceed 1 inch 


3 
Pl HPL 
6 2 2 eS Se (4. 4) 
eet ana 


Additionally, geometric limits are imposed on the Leam size. 


0.5 <B < 5.0 (4.5) 
10 See eo (4. 6) 
H/b < 10. (4.7) 


Now we can input this problem to MSCOP. 
Tho te NDY oN tOC ot ei 
UO 1T008data 2,42 


initial staremiGg pornmes 


Evaluation of objective 


00400 fn £ = tl*x(1) 4x (2) 
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Eraliiadtion Of COnStLraints 

oo 00 tl = 209. 

3823) Le 2 43588) 

o> 0 3 ae i = 1 then fn_g = Lee 
Piece: = 2eehem fnig = 3.%bp/( 10000. *2.*5b*h)—1. 
moo02 if i = 3 then £ng = ee ee eee he) 
00503 if i=4 then fn_g = h/E-10. 

TAGE 2 


The Solution of a Uniform Cantilevered Bean 


objective ; 6664.0 
design variable ; 


Meat) 
X (2) 


1.852 
Pe 


constraint ; 


g{1) = 0.000902 


g(2) = -0.9549 
C4 3a 0 0 109 
g(4) = -0.0286 


As a result of this froblem are in Table 4.1. 
2. Variable Cantilevered Bean 


— <a eece ee I el = —_ —— 


The cantilevered beam shown in Figure 4.2 is to be 


designed for minimum material volume. The design variables 
are the width b and height h at each of 5 segments. we 
wish ae design the bean Subject to Pe IS on 


stress(calculated at left end of each segment), deflection 
under the load, and the geometric requirement that the 


height of any segment does not exceed 20 times the width. 


Bo 





P= $0,000 N 
i £ =200GPa 
A, L=s00em 
6 =14.000N ems 
a ¥F= Sem 
—| b, |-— 
Cross section 
rc ee a er mh a a 
Figure 4.2 Design of a Variable Cantilevered Bean. 


The deflection y at the right end of segment i is 


calculated by the following recursion formulas 








Watartiy a4tl (4. 8) 
oO O 
Pei 1 
= = L + ee = 1h + (4.9) 
2 El, 2 31 i Tint . 
a 
2 
Dae 2 te 
= i: > Lt . + 1 (4.10) 
irae j= 3 7-7 en ; 
a 
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where the deflection y is defined as positive downward, y' 
is the derivative of y with respect to the X, and 1; is the 
length of of segment i. Young's modalus E is the same for 


all segments, and the moment of inertia for segment i is 


I = —— (4.11) 


The Eending moment at the left end of segment 1 1s calcu- 


fated as 


al 
=e [u+a,- | (4.12) 
1 jJ=1 21 


and the corresponding maximum bending stress is 
at 
¢ = —— (Cia) 
a 


The design task is now defined as 





N 
Minimize : = eee hee 9 (4. 14) 
a ear. 8 
Subject to : Cas 13) 
O-: 
= oat 1 < 0 aL = We terrereoiN (4. 16) 
oO — 
YN 
—— - 1<¢ 0 (4.17) 
yi 
he ae OD? <9 0 ep gore oN (tao) 
a i- 
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SO h > 5.0 i = 41 (4.19) 


where Gis the allowable bending stress and y is the allo- 
wable displacement. This is a design problem in 190 vari- 
ables. There are 6 ncenlinear constraints defined by Eqa(4. 16) 
and Eq(4.17), and 5 linear constraints defined by £q(40 lem 
and 10 side constraints on the design variables defined by 
Bo (ues 19 ja. 


Now we can input this problem to “SCOP. 
Input NDV, NIOC, IPRT 
OO TOC datas 0g sla 


Initial starting poernts 


0021 02davas ooo. iar nO 
00220 Gavtar >i. oO 
ODZS 0° Gata-sa., i 1.0 
OO2ZU0 data. s.5liesno 
0025 0 tdabay >. 5 1.5. no 
00260 data 4).,5.,n0 
QOZ70 data 402-5. ane 
H0Z280). data 10.55... bo 
UO290 “data, 10..5.).nc 
J0300 data 80. , 52 ino 


00490 fn £ 00 * ( x{1)*x (6) + x (2) x) + x (3) xg 
x(4) *x(9) + ) *x (10) ) 


( 


Evaluatiom of constraints. 


Ge ao bn 1%" it (10 19 bie b (10 
im b a JS af a We ) 

00500 pcb = 39 0b6. ; oe : 

00501 be = .et5 

00502 ti 5= 589: 

S050. 3 Sh = 74000. 

GOODC dy bbe 5 

OO 50s ‘Ste ——40'. 

O05 06 Ler am = 1 toss 

00507 bm(m) = pcek*(tl+sl-m*s1) 

00508 next n 

OO50S for m=] 1 to. 5 

00510 = km = mt5 

C0812 + (Shai (mye ba tm) *x (km) 7 (2. *bi (m) ) 

Sigqgi(m) = br xX (kn om 

OCS 14 next ; 

C0545 zo —— 0. 

00515 yrzo = 0, 

00516. for mea= Ween 


216 


DoOOOOQOCOOGDT Ooo OOVWOoO 
OOBOVWOOODOVOVOO DO VODO 
Or NnNnNnnnDd HOI W1N1 1 10 0 II 
WJ WLILWNI 32D OM~JOUINI NIN 3 
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b(m) = Guroit (tlL+si/2.-m*si)) 7/7 (be*mi (m))-ypzo 
ee z Beets 2 (tl-neGl 2. %2l/ 5.) eee 
yb(m) = dunmm/ (2.*be*bi (m)) typzo*sltyzo 
Ypzo = ype” 
ze = yb(m) 

next nm ; 

Een COiauearnt EUnCLION. 

iio — wieeneme ta g¢ = Sigi (1) /sigb-1. 

hel) = ceo ened = Si gi(2) /sigb-1. 

Poel mene nine g = Sidi (3) /sigb-1. 

1£ 1 = 4 then fn_g = sigi(4)/sigb-1. 

il —s oeemen fn G)= Sigi{d Aer oe 

1£ 1 = 6 then fn 7g = yi 15) 74 Dal. 

ie ene theg = x(6y—=20 ox (1 

if 1 = 8 then fn _g = x oo 2 

ir = Jeenen freq = x(S)-20.*x (3 

ff a = §Omehen in g = x ay ao ee t 

heel = eleene Meth 9 = 9 xX (10) >20.*x (5) 
TAS EE ars 


The Solution of a Variable Cantilevered Bean 


objective 


62133. 35 


design variables 


1) os 
%(2) = 
X(3) = 
a(S) = 
K (3) 
X¥(6) = 
X(7) = 
X(8) = 
(oe 
X (10) = 


2.994 
Pel O7 
223928 
222038 
de) Gil 
59.88 
Bio Ore 
Be D6 
44.14 
35.19 


a7, 


constraints 
G(1) = -0.090219 
G(2) = -9.00415 
G(3) = -0.00508 
G{4) = -0.00406 
G({5) = -0.0177 
G(6) = -0.4401 
Gq{7) = -0.0101 
Gleyw— —050231 
G{9) = 0.0000 


G(11)= -0.0278 


EB. WsdSPle ness 


ep = so0C0Os 
L 
E = 2Z00"Gra 
One + 14000 N/eom 
Xx tL = 100 cm 





—-—— ame eee ee 
eee = eee ee ee ee ee eee 


Figure 4.3 Design of a 5-Bar Truss. 


A simple truss with 5 members as shown in Figure 4.3 is 
-designed for the minimum volume. The design variables are 
the sectiocral areas of the members. The const railnts ame 
formed for the stresses of the members not to exceed the 
given allowable stress. The lower bound for each design 
variable is also considered. The stresses are oktained by 
the displacement methcd of the finite element analysis. The 


equation to be solved is given by 
K-u =P (4. 20) 


where Kis the stiffness matrix, uw iS the )dispiec na 


vector and P is the load vector as follows 
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u 0 
1 
Vv 0 
U = 1 p = Co o21) 
= u = 0 
2 
Vv aod 
2 
A As As 
L ee GL 0 © 
As A, As 0 eee 
¥2 2 “Q y¥2h Q 
K See (4 22) 
0 o As, Ae _ As 
Leer 2 
O — Aa ate oes 
2 x ez © 
From Eq. (4.20) the displacements are solved by 
-j 
eee owe: (ie) 
Having displacements at all nodes, we can calculate the 
stress for each element. 
E?Al 
1 
O;= E-g = (4. 24) 
i. 
al 
where 
2 2 
41 = j}(1 + uu) +Vv =e gt 
71 1 i 1 1 
Z 2 
41 = (eo = oye) + ( ue = uy = 1 
Zz Z 1 2 1 2 Zs 
2 2 
41 =f 1 +u) + Vv ee (4.25) 
3 3 2 2 3 


Si 


ae 2 2 
AI =/( 1 + yu)” + °C tee eee 
4 3 2 4 
2 2 
Al = K L+#u) + ( leven 
“5 a> 4 a 5 


The design problem is given by 


5 
minimize VS 2 eee (4. 2G 
i= 1° 
Subject to 
| 0: | | | 
So S10 eee (1.295 
a Oo - 
a 
BT DAO i =. lee (8.28) 
a = 


The MSCOF input for this problem is given as follows 
Inpuc ino)  NEOC. SIert 
COND OF dated. 57, 51,2 


TnVival sctartinGg eon 


IO2Z0 0° Gatak 3c co tno 
QO2Z02 “Gata 3.7.0) -n0 
00204 data 3.,.1,no0 
00206 (Odea 3.7.07 1e 
U020'8 “data 3.54.1, 00 


Evaluation of objective 


00409 fn £f = 109 * ( x(1) + x(2) + x(3) + Sar(2.)¥*x(4) + 
Scr (2.4) =o) 


Evaluation of constraints 


0500 “da tevwte) 
050) te = 2ret 
Veet = 1008 
0593 sigb = 14000. 
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TABLE III 
The Sclution of a 5-Bar Truss 


ase >< 


e 
e 


objective 


constcaint 


design variables 


-~1.9988 
~2.90030 
#0. 0030 
-). 1203 
=1207 97 


G (1) 


0. 1 


X (1) 


G (2) 


0.1 


X (2) 
X (3) 


G (3) 


3.514 


—= 
=e 


G (4) 


4.948 
Ont 


X (4) 


G (5) 


X (5) 
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Ve. SUMMARY AND CONCLUSION 

Numerical optimization is a powerful technicue for those 
confronted with practical engineering design problens. iets 
is also.a useful tool for obtaining reasonable solutions to 
the classical engineering design problems. Since many engi- 
neers are now using wicrocomputers for solving design crob- 
lems, the development of microcomputer software which can be 
easily used is needed. 

In this thesis, an algoritnm for constrained Outi 
tion problems is programmed in standard BASIC language 
(HBASIC version 2.0) senvan tS tesCcer The users can easily 
convert this to other microcomputers. 

MSCOE (MP CEOcenp iter Software ene Constrained 
Optimization Problems) employs the method of feasible direc- 
tions and specific mcdifications of a one-dimensional search 
£Or CENnSstraimed ‘oOpemnizacion. MSCOP has been validated bv 
tests on three constrained optimization problems. Tte 
performance is good and could be made better through refine- 
ment cf the algorithno. 

Since microcomputers are available with reasonable 
memory size and computational speed, their capabilities will 
continue to improve as more engineering software becomes 
available. MSCOP is considered to be a first step tcward 
more widespread use cf optimization techniques on microcca- 


puters. 
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APPENDIX A 
MSCOP PROGKAM LISTING 





cretion base 1 

Aa Scie Omeenneedeviol) ,ngcv (51), — ot} 
dit thta (21 ,we 159 51 : 

Geiss lh 2 | aeear sy peo eal) yen es 2 eGo) 
ait wec (aly. 1wrk (51) , wok 5 WIOK2 g 4 rer k 3451 

din sweirGowy, wok (51) ,lowb (21 ,uprb 21 on reece) 
rem input data 

GoStees00 00 


Pet eomout num. D er a design variables and constraints. 
read g Ne hove y pet 
data 
EOE ee z "4 to nay 
rem input initial value of design variables 
read x (i) 
ONG) — x (2,) 
if nigc = 0 then 160 


read lo$,up 
if loS = "not then lowb(i) = bnio else lowb(i) = 
yaeue {1 oF) . 
oe U "no’ then uprb(i) = pnup else uprk(i) = 


value (ups) 


Xt 
sea aa Sy Oe ager 
data jee, 1 .0, DRG) 
rem evalite the co vert ane tion 
obj.= a2 _f (x) 


ren objective function 
def fn f£ (x ) 

eee ek xX ( 1) * xX (2) 
fnend 


rep evaluate the constraints - 
fOrESi== tb to nigc 

ay (1) = fn Poi, ets) 
nex 
rem constraint functions 
see fn ee - 2} 


be = 49, e+6 
PP = 10000. 
een fn_g = eA 
J 2009 eae ac i 
eee le rng — (s.'*bp) 7 (20000.~x (1) *x (2) ) -1 
ties oethen £neg = (4. *bp*tl**3) 7 
. ae exe ie Se Si | 
Toot nen  f£neg = x (2) /7(10.*x (1))-1. 
fnend _. 
Peniitetal GOuntIng number input 
ical =_1 


ete d er 5 then Stop 

BeCimeall the OfPtimization code. 

gosub 2000 

rem print results. 

rem 

rem re-counting number input. 

ical =_ical+1 

if ical = 3 then 850 

rem 10% reduce the design variables. 
for 1 = 1 to ndv 
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LON NRNNNNNNNOCOOOOOCDO0009D 


ODWOOANUIE WIND 29 OOO NID N EWN B@OWD ANNE WYN @aOWOONNMNIEWN BOW OWNAN LW pO 2OOOdDOOWO~AINU1I4S WN = 
OODO DODO OS OSS SO SOOD DOOD OOOO MO VOD DO BO VSO DOSS OOSOOSODOOOOD OD DOOLWN  CDOCVVDIOOVOONOO 


ONnNSEeSSE SLES eS oS WwWwIOWOIWWIWIWIONIN NINA DARONIAA) 323 3) 25) 2 a OOOO COO DODOCDOOOODOOMO Bow 


WMDNINNDNNNNDNNNNNNNIVMNN DIN NNINNNNNDNNPINMNNNNNNNNDNNNN 


ee = 05 ox (a) 
x0 (= xn) 
mex G4 
Goto 7207 
Tem 107 increase design variables. 
for 1 = 1 to ney 
M (2 )e =. 16, 1 ee) 
MOS) Tee ee) 
next i 
goto 7/20 
rem calculate the obj. censtraint fen. 
obj = £n_f (x) 
LOL 19 ao 
fe (1)ic= Envoyer 
next i 
1trgq = 


LELd = teary 

rem calculate the number of active and violate 
COnStratn esr. 

gosub 3500 ; 

rem calculate the gradient of objective and 
active or violated constraints. 

gosub 3800 

Li fave = Urthen, 2170 

gosub 3900 i 

rem calculate the push-off factors 

gosub 4000 . 

ren Me eee a Bis te 

rem normalized the df (i) 

gosub 4100. 

rem normalized the DG(1i) 


gosub 4200 
if nvc_> 0 then gosub 4400 else gosub 4600 
rem calaulate the usable-feasible direction s(i) 
gosub 5000 
GOtou22 0 
rem normalize the df£f(i) 
for 1o=. 1 torndy 
Se = a ee 
Next 7 ; . 
rem normalize the s(1) 
gosub 5700 


rem one-dimensicnal search 
if nvc = 0 then gosub 6000 
Len vupdate x fOr ain 
gosub 7000 


else gosub 9000 


-accx and delf <= dabf then 2470 
bd = 9 De eee . 
tee ae then print ‘check the prop lem 


forygs. — | torvnaees |. 
gov ta = fn_g(x,1) 
next i 


LE UOrt.=~2 thenezeoo 
gosupe 9.20 

Touro 72 : 

Lem Print, final resulee 

Print **#A4* fanal results. eee 
gosub 9200 

return. . . ; ; 

rem initialize the integer working array 
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MUNIN NIE PSHE HPS LESS WWW WWW yw WWW WWW WII DN DDD poh phd phd —2 2 2 ee ot er OOO OOOVOOOD 


Wh wOQ ODOYVIANQAUIWNIW NN 22 OOM ODO YOUN OAHUNUIUITUIN — 3 OOO YI OAQHUIIN 3 2 OOD DNDN 2 OO LINHA a 2 OO 
ODOLSN ODUODNDNIOOUVWNG AO' NWO YH AaADUINUNWWDONODNBDNONDNIVD IO NID NWDOWNNIO NDA NWNIBDIOWIIODNYDVIOWN1IGW1 


font = ) too nacgn 
iwrk(1) = 0 
next 1 
return 
rem initialize the integer working array 
GO eee coun. om 
a) = 0 
nex 
EeCuUrn a ; 
rem 1eypak ese) Wak vAcs 1eiel= one-dimension working array 
feu 4 = to nigm 
wok1(1i) = 0. 


i eooGeL 
return ; 
rem initialize the one-dimension working array 
fore) = Ietko nign 

SEK = O. 
next 1 
return ; : 
rem initialize the one-dimersion working array 
SOR = | tO ace 

Wek stl) = TgGy (1) 
next 1 
return . 
rem initialize the two-dimension working array 
£Or 71 =.1 to nagm 

£Om e) = i terwndvm 

woky (1,3) = O. 

next 
next il 
return 


rem anutiailaZze the derivative of objective DF (i) 
fEOmwt —mei tO ndvm 


Gia(iyee=| 0. 
next 1 
return 
rem initialize Ens leew Dil) op aajagc (2) 
£OL i= | fe ndv 
1) = Q. 
Pi} = 0: | 
Or J =,1 tc nigm 
ayia, ee. 
Peace. 3 
Rescue 2 : 
One = 1) FO Dicm 
14) = 0. 
next Jj 
return, . . 
rem initialize the derivative of constraints D€(i,j) 
Ort — | tO nigm 
FOEgjme= | to ndvm 
agin, 3) = 10. 
epee F 
next 1 
return 
rem initialize the b(i,j) 
FOE = 21 Commi qn 
fGEagie— 1 tO 210m 
b(i,j) = 0. 
next j 
next 1 
return 
rem Calculate the number of active and violate 
constraints. 
gosub 3090 
gosub 3700 
hac = Q 
hve = 0 
Porm. tO Nigc 
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~ADOO LS BWWPNN 9.2 OO OU SWNwaODnHDINIES SWdntle saOONUNEWWDIN @ wrODWUIE WN 2D OOADNDUIE WN 2 OD~) HONE 
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1 Cie: >= vec then 3580 
aLeCey < SSG then 3590 
nac = nact’1 
goto 3550 
hve = nvct+? 
next 1 


have = nactnv 
it Mane 5 0 then a7 2c 


ii = 
Soa 
E.2 = | te nige 
1£ gcev(i) >= vee them 2720 
1f gQCcV(i) < ace then 7a 
iwrk(nvc+ii) = 
wrk tae eee = gcv(i) 
12 =" 2a 
JOLro, 37 50 : 
LW ete 
wok1(9) = gev (i) 
25, = Jo! 
next 2 
return 


rem calculate the gradient of f (x) 
osub 3300 

Or 2 =I to: may 
dxi = £fdm*abs (x (i) 
dxi <= mfds then dxi = mfds 
i = x 2 +AxX1 


= Ager obs) 7axi 
x 


~~ OO b- 
Kh | 


1 
idk.) Sets Veiga dxi = mfds 
x (1) a oe) +dxi 
J = eG Mave 

kv = wrk (3) 
dcon = Fng oe 
one ( Sone WOK 1 (3) }aet 

next 

Kea) =o xOn) 
next 1 
return 
ren ec aes eae puUSh-OLEVraecter 
for i= to 

thta (4) = t hto* (1. —WrK 1 (1) Jace) 4.2 

see thta (i) > thtm then tht aauiee= sce 
next 1 


return 
rem normalize the DF(1i) 
ogee a0. 3200 
om i = 1 to nd 
fSen = fsqtdt (i) **2 
Pe af f 
Sg = sgqr(fs 
eee 2 cP then fsq = zro 
fous to 
oe = 13 »/£8a) *a€ (4) 
next 1 
LectuEn 
rem normalize the DG(i) 
osub 3250 
Ob A =") tO eave 
gsg = 0. 


YE 


OWANIA AAAI AUIS ek ee eer ee PE ee EERE reer EePeeEe eS eereEES eS eee eEerL ere 
34 3 2 2 3 BOO ODO O0 000 ~~ JJ J JIS ANNA HANAAAAAAHMAINAIMNMANNANININ SSeS REE EEE EENNNNNMNNNN EY 
SHANE WN BOW DYANE AOD ODNNUIEWNRPOUOAYDNE WN 2ODUOOONDNHDNEWN 2OVDOYNHNIEWN2OOWNUIE SWWWNN 
DOD DOO O DO OVO OOOO OVUINIAD DOOYOO OOD 0@ODYOO OO DVUIO DNOO OVO OO BOD OOO VD DOOOMOMNIONOUVISO UND tow 


hor) elec 1a 
S 


e qi gsq+dg (i, a2 
next j 
934 = 22S a et - 
gsq en gsq = zro 
ae J oenday = 
weky (i, 55 = (1./g9s8q) *dg (1,3) 
next Jj 
next i 
return 
rem exist the violate constraints 
osub 335 
Of. — 1 ro Nave 
Cee eee iC ven 
a 1,]) ir wrky (1, J) 
next Z ; 
an iy, eev oo) =e td (2) 
Meat 2 
f Olea = yr tO ndv. 
p (2) = -wrk3 (2) 
nex 


1 
p (ndv+1) vhid 
or i= 7 to nave 


br 1 tc ndv+1 
a (1,)) *p{j) 
YY+XX 


ms 
he Ps, 
1 


J 
Ai) ees Vo) VY 


osub 3 
Om 3 = 1 Mg havc 
44° hema  - 
a(i, 5) sew ay (79) 
nese J 
a(eomayv+t) = tChta(2) 
mete a 
for j = 1 to ndv . 
a(nave+ 1,3) = wxk3 (3) 
next. 9 
ee = 1. 
P ndv+1) 1 
oie 2 
SG 


Oe 


to navct+1 
a(i, Bea) pny 1) 
Gj = 


(-1.)# 
next 1 
ndk = navc+t1 
return : : . 
rem calculate the usable-feasible direction 


f a fe et (k, J) 
= atl Wr . 
f = paste sete 


= (-1.)#ff 
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next 1 
iter = 0 
nmax = 5*ndb 
rem begin iteration 
l1ter = iter+’ 
ChuxX = OF 
ichk. = 0 
fOr 12 = Vto nae 
C1 = c(i) 
Dai dese ee} 
if bai — 0s sthemeoeue 
if ci > 0. then 240 
Cb = C47 pi 
Pi Chige ebnx then 5349 
ichk = 1 - 
ckmx = cb 
Ne <<Go . 
if cbmx < ZrO OF iter > Dax then.) 
if ichk = .0 them ao a0 
Vis] AWEK(1ehk) aa 
ee 34 = 0 then ee ee = ichk else iwrok (iL chk) = 0 
i= b(ichk,ichk),.= 0. then D(izenkP7 sent) 
DD =. 1.7 2 Gerwarenk) 
1f bb = QO. then bb = zro 
LOL fi s= 7) found 
b(ichk,1) = Db*b (teh 
next 1 
€{1¢Chnk) y— cbhnx 
LOG. 1-=) 1 ton nab 
ifs i= Ten ethen o550 
DDie= oo (i, a-eiks) 
b(a,ilenk) = 0.2 
£OU j =.Ameo ndb 
dei 2 = Lenk then 55209 ; 
(i,j) = b(i,j)—bbi*b (ichk, 3) 
next 
ene) = c ti) -bbi*ebmx 
next i 
GOtOmo 220 
ner = 0 
FOr 4. =) | toende 
u(i) = 9. 
A] eke 
aia 2 > O-then (1) = ¢ (3) 
next 1 
POmeto => is GOsn dr 
fin = O. 
Hor 1 = nd 
PE p= Et tw ky (a) 
next j ; 
i) 2 ere 
S$(i) = y(i 
Hee 
return ; ; 
rer normalized the s(i) 
ssq = Q. 
On. 1 = 1 to nev 
SQ = SSqtsta)=-2 
nee 2 


SSsg = _sgr 
aeeEe pee ie a hen fSip — Zino 


fOr w= 1 toewndy 
S(i).=°41. 7seq) =a) 

next 1 

return 


rem one-dimenSicnal search for initial feasible pcint. 
rem calculate for slope of £({(x) 

fslp = 0. 

Forel = “Veto endy 


48 


NOM ANOWONOSNONONO NON ANOUVION MNO VAN NBII 


HADDAD DARADAHAARNHAAARAADAHANANHAAAANANAAAAHARAA DAAHAAHAAANRAAAHANAAAAHAAANANAANAANANAARDAAHAAAAAANA AO 


WIILW WUILILINDODNIDODI DODD AODNONINIDIDODONIDINIAQDINIDININAND on a 3 a 9 2 OS Ss Or BOO ODOOODOOOCOOCOOCOOOCOOO 
RONIDN SRM ODOWOUODOOAYDINNUINUIE FS WWWDNNAIOO YOWOODO~YDWINONUIUNE SWWpohddd 3 2OOW OO D~)~JAION OHA & Woh) 


HANIa QD NOVO NOWNOUIO (ADVUIN OW10 YN AOUIODUIO 


fo) Pero tar (1) *s (1) 


eee 1 « e e ° 
rem idenfy the initial point. 
alow = 0.. 
flow = obj 
POmel— sal bOmlle 

WE KING), = “-gev (2) 
next 2 


rem find alist ; the ist mid-point. 
ioneeamecwj oy — = (0) then Hee = Zro 
alst = Se Honan ( Sl p) 
EOE 2 = Siero ndy 
if s(1i) = 0. then s(1) = zZro 
walp = ee eek ey) ) 
lf Waip > alst then 6095 
_alst = walp 
He xt. 
rem update x for alst. 
alph = als 
gosub 7000 
gosub 7100 
rem calculate the flst and wrk1(i) 
fist “= one © (xX). 


EOGtow = 9) to nice 
wok1(i) = f£n_g (X,1) 
next 1 wie 
rem check the feasibility. 
ncev1_= 0 
Comets = ly eo. 2ge 
Pe werml(tha< vec then 6170 
nevi = nevitt 
next 1 
if ncevl = 90 then 6200 
ast = 0.52a 1st 
Goro 6 105 


MReimalnGeea 2nd =. thie 271 mid-point. 
Boies poOrines —GuUadratic Lit 1nterpolation 
Pore nine m £(aipa) . 


aO = flow 

al = f£slp 

a2 = (f£1st-al*¥aist-a0) /(alst** 2) 
fF a2 <= 0. then a2 = zro 


ae Gm a aie as) 


rem 2-points linear interpolation for g(alpa)=9. 
POEel. = 1 to nige 
ad = wrkl(i) 
medicte— 0. then alst = ZrO 
al = eo eeu a tst 
1£f al <= O. then al = zro 
walp = -a0/al 
1f walp <= 0. then walp = 1000. 
ae Ee >= a2nd then 6265 
a2nd = walp 
exe 
rem update x for a2nd. 
aden — a 2nd 
gosub 7000 
gosub 7100 
rem calculate f2nd and wrk2(i) 
Pde seit (%).. 
fOtet = abo N1Loc 
Wie) = ing (x, 1) 
next i 


rem find final ee eaee, by using 
SOO Mesmevatiatic fat, 

£1 = flow 

pole - alow 


Ist 
alp2 = alist 
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£35 =f 2nd 
ales = a2nid 
Een os Be50 quadratic fit interpolatiog. 


690 
if a2 = 0. then a = Zro 
= wrkl 
i205 We 1} 
_. gosub 6630 
li alps asd eee 6380 
ren update x £Cr auc 
alph =_a3rd 
rem calculate the fupr and wrku(i) 
EVEL ceee a EO 


ayes) 


Pf alJed <0 Eee a3rd = 100ge 
wrk2 
aged Saeed 
gosub 7900 


1 to ute 
gosub 6600 
next 1 
gosub 7100 


Leto Hence 
eo) =i Ng (xe) 


next 
rer find 4th new point. 
ae = f2nd 
= 2G a 
eed =atlst 
aleZ <= "ana 
alp3 =a3r 
rep arcane quadratic fit. 
gosu 60 
af az = OQ. then a2y- 256 
aupr = goa) 
£OEG1 =] i vosmig 
£1 = wae 1a 
£2-= Wek 2 (i 
£3 = WEKO 
al l= a list 
alp2) =2a2nd 
alps = asrd 
gosub 6600 
gosub 6630 
1f alps > aurr then 6540 
aupr = alps 
next 1 


aoe update x for aupr 


h 
a oean ~7900- 
gosub 7100 
rem evaluate furpr and wrku(1i) 
Ue eel (0) 
£Or i= 1) tO unwtee . 
ona (i) a= oe Ge) 
next 1 
rer find Beye oe alpa for not violating constraints. 
gosub 143 
return 
rem 3- oints Wadratac fies 
ey aa = a: cr alp2 = alpam@orn ale!) — ake 


then 
a2 = ({E3- £1) 7 (al Se a5 Ry ee 

2 2 £17 aiBs =a le as ane ce 
al = (£2-£9) 7 (alp2-atae ioitalne 
a0 = 6 Pao peat bee 
return 


rem zero of polynomial for g(alpa) 


aC 


09.00 CO 00 COCO 0D CO 00 00 00 09 00 0900 CO 00.00 09 00 COO SJ SJ MII SS I SI SS SEDDON AHNANAHA AN DANDANAANAANANANAHAADNAAADAARAARAAAN AO 
NON 2 3 3 SS 2 S 9 3 3 OD 9 DO OOOO oA @ 2 29NO0 OYMO\0O 0 WO 000 C0000 OO ~] ~} II SI J SIN HAKAHANDAAARARAHAATDNA 


=2OW DANN EWN 2 DW OWNING WN Rp ONES WN 290 SUN 2 OO WOT DUIUTUIN = 9 OW WOON 2 2 2OOWWODOYVAINNNUIE FE Ew 
DO OODDYODOYDOO DDOOO OOO DO VOOOO9O OOOO OOO OO OO Vi NDOONU1IDI1O OO MN DIO NOUNDNIDNOADUINIO V1 


dd = al**2-4.*aZ*a0 
Pema < 9. then 6715 
Le 
a 


Eo then 6710 
apie 


alrc 


ONNNHOANN Noow 
NO 


ON Nh 
© ONnowvoil 


return 
ren Dee a oaE 3b4 alpx 
delf ob fy 
Ger abs de £796] 
Aae Oj+diGay72. 


OF 


fw 

a: 

oO i 

LO. 
Wout a 


ausex = (alpxt+walp) 7/2. 
dabft = Wegeea be to0g) 
return 
Eemeupdat Ec x (1) 
for i= 1 ndv ; 
x (i) = Ocoee saan 
next 1 
ew) ; 
rem check the side-constraints. 
HOnei-= | to ndy 
Me (ee Bevo) then x (i) = Lowb (1) 
PEC (1 etehen x2)" Wor D1) 


rem estimate the alpa 

fstr = flow 

alpa = alow 

nvcl = 0 

BOr te — 9 | tO eae 
ete wr k1 (1) vec then 8079 
nvctl yon 


then 8120 
ae then 8120 


Ss 
st 


a 
rh 
th 
aod 
Nn . 
ct 
luvy 
thm FhO 


— — (/) 


EOE pi: peo nage 
at we k2 (2) < vec then 816) 
nvel = nvcl+1 


next 1 
O then 8210 


env cl 
ie oe ona fstr then 8210 
a2ndi 


fond 


nuvy 


Sal 
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NN NDIDNNNDN OO DOOD VOODOO DOOOO O09 OOO C8OOO DOWOOWOOODVWOAO OOOO OC LWWWWWWWWWWDNPNPNt Unb 


WIDIN) 39 OO CO YJ YAY ASIAN AANNNOINIONVE & EE EW WWWNIINNN YAMS Sw OOOO OCDOWODANNDAUNIL WIN 2OUWO OWN NE Wp 
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EOL Bel 4 le On) ae 
lf wrk3 (1) < vec theneeveg 
Hvel = nye 

eee 

if nvel > 0 ther 3300 

if £3rd > Estr then 3300 
aoe = a3rd 
fstr = f3rd 

nvcl = 0 . 

{Orsi =) (eto sn de 
if wrku(i) < vec then 8340 
nvel = nvcoltt 


next 1 
if pycl > 0 threneoso 
if £Upr > EStr thengee 6 
ae = aupr 
2 ESte, > eons 
aa = alpa 
re urn ° ° ° . ° 
rem one-dimensional search for initial 
_. snieasibte former. 
ae ee 
qcvm = wrk1(1) 
LOL 1. = 1 tornave 
LE wrk1(i) <= gcvm then 9014 
ae ee 
ee: = wrk (1) 
next 1 


rem calculate the slope of badly violation. 
slp. = Q. 
OL 3) =" bo andy). aa ; 
ae = 9st p49 (2a eer) 
next i 
rem calculate the alph. 
a3 85 aoeu = 0. then gslp = zro 


al -gcvn/qsl 

on fede X coe cents 
gosub 7000 

gosub 7100 


rem evalute the objective and constraint. 
Ob} = inf (y).-. 
for 2 = | tounagc 

as ae fn_g(x,1) 
next i 
rem calculate the NVC. 
gosub 3500 
1f nvc_ = 0 then return 
rem update initial value. 


POE 1e> 1. tO dy 
OC) = x (2) 
Nex, a : ; 
rem calculate df(1),dgq(i,39) and@push-oft -acitom 
gosub 3800 
gosub 3900 
gosub 49000. ; am 
rem normalize the df(i),dg(1,j) 
gosub 4100 
gosub 4200 ; 
rem find the search direction. 
gosub 4400 
gosub 5000 
oto 9000 
rem Poa the results 
DOAINC "#KHKKKKHERE GQtQq KHKKKRHKEKKE! 
print '! ; 
print ‘The number of design Variaqnve-wee = "gndvy 
pean \ the humber of inequality constraints = "nage 
PEM 


eZ 


OO OWOWWM ODOWW'IDODWO OWODWUDOWMOUWODWOUWWOWUWWWWOWWWUWWWOWWUNUW 
OY ADKDDN AKDAKATDHAN NAINA IWWWWONNNNNNNNNNPN 
OD OOWNN NMNEWNUDOO DANN EWN OO 2DOWOWOOD-JINDNHMNUIE EW 
OD CXOOO DODVODYD COVCSOODDIOONOUNONDWNWNOMNOWMNOtNOv1 


print ‘The objective vaiue = ',05j 


print : 
mine on cecaldgn Varlables *****" 
Or 1.= 1 to ndv 
pea te SC ee), Sax (0) 
next 1 
print (! . ; 
Berne ‘the number of actave constraints = "j;nac 
rin 
brint ‘the humber of vionate constraints = "3;nve 
real tc 
eae loa CONS rLalnt value ++" 
reais 9g eee : 
Or 1. = Metornidc 
pea e OMG 25 eee ey (1) 
next 1 
return 
rem default numker 
Hike = 50 ! maximum iteration number 
fdm = .01 ' finite difference oss 
mfds = .001 ! maximum absolute finite difference step 
vec = .008 ! violated constraint criteria (thickness) 
acc = ~-,1 PVacetvenCousesarnts Criteria eee bee a 
Eh EO: = ae | pushs=ett factor Miltiplier (theta zero 
tiv ae Maxum value of push-off factor 
phid = 100000. ! weighting-factor used in direction 
when infeasibie 
aeecr = 200) ! absolute convergence criteria 
accx = 0.001 ! absolute convergence criteria. 
ZeeGe= «001 ! defined zero a 
espl = .005 _! used to een division by zero 
bnlo = -1.e+70 ! the value of low boundary 
Dip — tect, 0 'Sthe value of upper boundary | 
dalp = .01 ! step size of alpa in one-dimensional 
Sealene . 
abicum=. 0. | 1 step size for reduce objective 
alpx = .1 ! reduce the design variable factor 
ndvm = 21 ! tke number cf maximum design variable 
nigm = 51 ! the number of maximum inequality 
CenSpra wits 
me Clean 
end 
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